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invariant-mass correlations and the resonance 
production in a thermal model with single freeze-out. The predictions are 
confronted with the recent data from the STAR Collaboration. 



Detection of hadronic resonances at RHIC (see the contribution of Fa- 
chini and Refs. H Ej) provides important clues on the particle produc- 
tion mechanism and subsequent evolution of the system formed at mid- 
O I rapidity in relativistic heavy-ion collisions. In this talk we analyze the tt+tt" 

invariant-mass correlations in the framework of the single-freeze-out model 
of Ref. ^ . For more details the reader is referred to Ref. 0]. 

The thermal or statistical approach to heavy-ion physics has been very 
^ ■ instructive and successful in classifying and understanding the data (c/. 

^ , contributions of Calderon, Gazdzicki, Hama, Kisiel, Ster, and Becattini). 

We recall briefly the basic assumptions of our model: 

1. An approximation of a single freeze-out, Tchem — 2kin = is made. 
This has recently gathered experimental support, with the short life- 
time of the hadronic phase obtained from the HBT measurements 
of i?out/^side ^ 1 from the out-of-plane elongation of i?side (c/. 
contributions of Appelshauser and Lisa). 

2. A complete treatment of resonances, with all particles from the Particle 
Data Tables incorporated. 
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3. The shape of the freeze-out hypersurface is assumed in the spirit of 
the Buda-Lund model 5 . It possesses a Hubble- hke flow, with the 
four- velocity proportional to the coordinate, = x^/r. 

4. The model has altogether four parameters. The thermal ones, T = 
160 MeV and = 26 MeV at y/s^ = 200 GeV, are fixed by fit- 
ting the ratios of the particle abundances. The geometry parameters, 
namely the invariant time at freeze-out r, which controls the overall 
normalization, and transverse size Pmax are fixed by fitting the p± 
spectra. 

In order to analyze the tt^tt" invariant mass spectra measured in Ref. [21 
with the like-sign subtraction technique, we make the assumption that all 
the correlated pion pairs are obtained from the decays of neutral resonances: 
p, K^, (jj, 77, 77^, /o/cr, and /2. The phase-shift formula for the volume density 
of resonances with spin degeneracy g is given by the formula [HI 

dn _ r d?p d6^^{M) 1 

dM ~^ J (2^)3 TvdM f y/M^ \ 

exp I y I — ^ J- 

and was used in Ref. 0j. We note that the same formalism has been also 
applied by Pratt and Bauer T in a similar analysis. In some works the 
spectral function of the resonance is used as the weight, instead of the 
derivative of the phase shift. For narrow resonances this does not make a 
difference, since then d6{M)/dM 2^ 7v6(M — mji). For wide resonances, such 
as in the scalar-isoscalar channel, or for effects of tails, the difference between 
the correct formula and the one with the spectral function is significant. 

In order to get a feeling on the problem, we first carry on a warm-up 
calculation, where a static source is used. Fig. 1 shows the mid-rapidity 
invariant-mass spectra computed at two different values of the freeze-out 
temperature: T = 165 MeV and 110 MeV. We note a clear appearance 
of the included resonances, as well as a feature of a much steeper fall-off 
of the strength in Fig. 1(b), which simply reflects the lower value of the 
temperature. Thus, the tt+tt" invariant-mass spectrum provides yet another 
thermometer for the decoupling temperature. 

The collective flow has no effect on the invariant mass of a pair of parti- 
cles produced in a resonance decay, since the quantity is Lorentz-invariant. 
Nevertheless, it affects the measured results since the kinematic cuts in an 
obvious manner break this invariance. We have performed a full-fledged 
calculation in our model, including the flow, kinematic cuts, and decays of 
higher resonances. The results are shown in Ref. 0. One of the thrilling 
aspects of the analysis is whether the observed shift of the p peak can be 



3 




^ ■ M[GeV]' 

Fig. 1. Calculation of the tt+tt" invariant- mass correlation from decays of reso- 
nances in a thermal model with a static source, (a) The case of Ttherm = ^chem = 
165 MeV. (b) The case of Ttherm = 110 MeV. The numbers indicate the position 
of the p peak (in MeV) and the labels indicate the decay channel. 



explained. We have found that the thermal effects and the kinematic cuts 
are able to shift it down by about 30 MeV, leaving a few tens of MeV 
unexplained. For that reason, we have repeated the calculation with the 
position of the p peak moved down by 9%. In Fig. 2 we compare the model 
predictions to the STAR data. Our results were scaled by a common factor, 
which is equivalent to fitting the value of r for the peripheric collisions of 
the STAR experiment. Then they were filtered by the detector efficiency 
correction (we are grateful here to Patricia Fachini). We note that the 
model does a very good job in reproducing the gross features of the data. 
We also note that the full model, with feeding from higher resonances and 
with flow/cuts, which has T — 165 MeV, gives similar predictions to the 
naive model at T = 110 MeV rather that 165 MeV. This is yet another 
manifestation of the effect of '^cooling^^ induced by the resonances 8 . 

More details and further results and discussion concerning the abun- 
dances and the transverse-momentum spectra of various resonances can be 
found in Ref. 
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Fig. 2. Single-freeze-out model vs. the data of Ref. [2j for the tt+tt" invariant- mass 
correlations (from Ihttp : //www . star . bnl . gov/^pf achini/KrakowModel|). The 

model calculation includes the decays of higher resonances, the flow, the kinematic 
cuts, and the detector efliciency. 
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